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Abstract 

The factorial moments (FM) of multiplicity distribution are used to 
study the deconfinement phase transition in heavy-ion collisions. The 
relation between FM and the partition number, M, results positive in- 
termittency exponents, 4> q . According to the signatures suggested from 
certain statistical models, the two-dimensional results of the dependence 
of <j)q/4>2, anomalous fractal, d q /d,2 and Renyi dimensions, 1Z q /lZ, and the 
normalized exponents, f g , on the orders of FM evidently supply evidence 
for the quark-hadron phase transition in Pb+Pb collisions at 158 AGeV. 

PACS: 25.75.Gz Particle correlations, 05. 40. -a Fluctuation phenom- 
ena and random processes, 64. 60. A Fractal, and percolation studies of 
phase transitions, 25.75.-q Relativistic heavy-ion collisions, 12.38.Mh Quark- 
gluon plasma 

1 Introduction 

The non-statistical fluctuations in the final state of particle production are 
suggested as an experimental tool to diagnostically confirm the elusive quark- 
hadron phase transition [El 0, B |J. Ever since the observations of spike-events 
first observed in the cosmic ray interactions Q and later re-produced in the 
laboratory || and since the pioneer works of Bialas and Peschanski Q, the 
intermittent behavior has attracted a lot of attention. Therefore, it has been 
examined in different interacting systems. If the quark-hadron phase transi- 
tion is to be produced in heavy-ion collisions, its critical aftereffects (like soft 
and hard collisions, clustering and resonance decay, showering processes, criti- 
cal exponents, etc.) are expected to survive even until the freeze-out Q. They 
reflect themselves in form of non-statistical (e.g. dynamical) fluctuations, power- 
scaling behavior, self-similar branching, mono-fractal density fluctuations, etc. 
In this letter, we study FM of the multiplicity distribution through successive 
partitions in one- and in two-dimensions (pseudo-rapidity, ry, and/or azimuthal 
angle, 4>). Using these investigations, we shall try to access the deconfinement 



phase transition 0] . For this destination, we will utilize the attitude of intermit- 
tency exponents, anomalous fractal, and Renyi dimensions on the orders of FM. 
QGP signatures principally suggested according to certain statistical models will 
be read off by means of the relation between intermittency and multi-fractality. 

The data sample used for this work is retrieved from some of the Pb-chambers 
irradiated at CERN-SPS during 1996 for the EMU01 collaboration. Specifically, 
the data sample employed here is that we have completely analyzed by our au- 
tomatic measuring system MIRACLE Lab Hat Marburg university. For details 
about the emulsion irradiation, we refer to (lfj . Due to their high resolution and 
geometrical acceptance, the nuclear emulsions are effectively capable to measure 
the charged particles, the angular distribution, and the density fluctuations of 
high multiplicity collisions. The scanning efficiency in emulsion chambers is 
~ 0.75 ± 0.05 |pl| . The sensitivity for singly charged particles is as good as 
30 grains per 100 /im. It is therefore, possible to locate the track positions 
to an accuracy < 2 /jm. Depending on the topology of the microscope's field 
of view, the produced particles with space angles, 9 < 30° (pseudorapidity, 
7/ = — lntan(0/2) > 1.32) can be acquired. 

The produced particles are expected to be mixed with contamination of 
electron-pairs from Dalitz decays and undetected 7-conversions. The possible 
overestimation of particle density has been determined as ~ 2% As a reason 
of the reconstruction algorithm applied for MIRACLE Lab, the tracks of these 
electrons are completely disregarded or percolated. The efficiency of MIRACLE 
Lab is estimated as ~ 96%. From manual/automatic comparison we noticed 
that the automatic reconstruction underestimates the multiplicity, especially 
the particles with relatively wide angles (77 < 2). Besides these missing mea- 
surements, the frequent scattering, unresolved close-pairs, nearest-neighbor par- 
ticles, and pair production represent an additional source of this 4%-discrepancy 
H . Many of the extra tracks in the manual measurement are close pairs sepa- 
rated by < 1 /im. The two track resolution, close neighbors, and the split track 
recognition are discussed in |J. In the automatic measuring facility, the lake 
of close-pairs within < 1 /im can be referred to the limited instrumental pair 
resolution. Also if the particle trajectory is spilt into two tracks, which can 
result a huge intermittency signal, both of them are rejected, since their own 
vertex definitely will not be coincident with the common one §. Then we can 
summarize that the data sample grabbed by MIRACLE Lab is evidently liber- 
ated from the double counting, measuring bias, and from any contaminations 
from particle decays or secondary interactions. We could therefore, consciously 
suggest to disregard the effects of Dalitz decays and 7-conversions on FM M 



(see Sect. 2.2.1 below). 

The accuracy of measured rj and <f> is practically depending on determination 
of the event axis and on the emulsion plates, into which the corresponding track 
penetrates until it leaves the field of view. The x- and y-coordinates are ideally 
to be measured with respect to some reference points whose positions are well- 
determined. The location of track in such transverse plane has a statistical 
uncertainty of ~ 0.3 /im. On the average, the tracks leave the field of view at 
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distance ~ 75 /j,m from the event exis. Then the resulting uncertainty in r\ is 
~ 0.004. Due to small air gaps, non-uniformity of the spacer between the plates, 
variations in emulsion thickness, the ^-coordinates are also uncertain. For most 
of tracks, the last measured emulsion plate is located 2.5 cm downstream. The 
uncertainty in <\> is then ~ 0.136 mrad |L2|, 




Figure 1: Pseudorapidity distribution of central Pb+Pb events 
measured by MIRACLE Lab (solid line). The dotted line rep- 
resents the distribution of FRITIOF, meanwhile the dashed line 
gives the Gaussian fit. The FRITIOF sample is simulated with 
zero impact parameter and default configurations. 



Fig. |T] shows the pseudorapidity distribution of Pb+Pb central events with 
multiplicity > 1200. In order to compare the experimental data to the ex- 
pectations based on incoherent models, we have simulated a sample of Pb+Pb 
collisions using FRITIOF 7.02 Monte Carlo code with very small impact param- 
eters. FRTITIOF code has been run in its default configuration. The dotted 
Line shows the pseudorapidity distribution of such sample. For 77 grater than 
~ 1.8, the two distributions are in good agreement. There is a small flattening 
in the central peak of experimental distribution. The flattening is expected, if 
QGP had been produced [^3| . 

If a phase-space of width A is split into M equal bins of size, 5 = A/M, the 
scaled factorial moments of multiplicity distribution are defined § as: 

Fq (M) = M ^ f; (i) 

n m is the multiplicity in m-th bin and n is the average multiplicity in the whole 
A-window. According to the self-similar density fluctuations Q (e.g. QCD 
parton cascading), the successive partitions lead to the following power-law 

§@: 

F q (M) oc M*«. (2) 
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The term "intermittency" used for the turbulence in the theory of chaos to 
describe the development of a hydro-dynamical system from stable to chaotic 
state can analogously be used in the particle production jl5|, The exponents, 
4> q , called "intermittency exponents", can be determined from the asymptotic 
behavior. Their behavior elevates the interpretation of the intermittency as af- 
tereffects of quark- hadron phase transition [^|, [l7j Based on 2D Ising 
model [^0|, [2l], ^2j, the same interpretation is obtained. In addition to these sta- 
tistical results, it is proposed that the intermittency parameters possibly contain 
signatures for the phase transition in heavy- ion collisions fi E|. In the contrast, 
assuming that the self-similarity dominates the particle production, it was diffi- 
cult to concretely suggest their intermittency as signature of QGP ||, |l4j| . This 
present work distinguishes between all these scenarios. 

4> q are related to the "anomalous fractal dimensions" though the following 
relation p^| : 

d q = <t> q ■ (q - I)" 1 . (3) 

At the critical point of Ising model |22| , d q can be given in terms of critical 
exponents. Also at this point, d q are independent on the orders, q. Generally, 
the q-dependency of d q interdepends on the nature of processes taking place 
in the interacting system. Therefore, d q are liable to the branching processes 
that precipitate an intermittent behavior. For this reason, d q can be successful 
candidates to give a further signature for the QGP formation 

The fractal Renyi dimensions, which ordinarily, are used to measure the 
randomization in particle production, are, in turn, depending on d q ||, |24|, |25| , 

n q = n-(i-d q ), (4) 

The constant, 1Z, represents the topological Renyi dimension. For multi- fractal 
processes, d q are linearly depending on q. For mono-fractal density fluctuations, 
d q are constant. 

Assuming that the intermittency and the mono-fractal density fluctuations 
survive the further phases until the freeze-out ||, different proposals are sug- 
gested to indicate the deconfmcmcnt phase transition via the interplay between 
intermittency and mono-fractality. From statistical models, it is expected that 
at T c the intermittency are related to only one particular combination of the 
critical exponents. In Eq. ||, <p q have been squeezed out in d q and then in !Z q 
(Eq. |]). Hence, the parameter governing the intermittency at a thermal phase 
transition can be represented by constant lZ q . Ordinarily, it is supposed that 
the mono-fractality measures the random sets of intermittent clusters of the 
confined phase inside the deconfined one ]26|, |2lJ . In a second-order QCD phase 
transition, it is expected that intermittent fluctuations take place and therefore, 
are characterized by unique anomalous fractal dimension. In 2D Ising model, 
the same behavior with constant d q has been proven |2^, ^6) . The validity 
of Eq. H is proposed as a signature of QGP 11), i.e. the validity of d q jdi = 1 
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according to Eq. ||. The expectations given in this section represent the model 
used in this letter to study the quark-hadron phase transition. 



2 Analysis and Results 
2.1 Scaled Factorial Moments 
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Figure 2: The q-order FM are given as functions of the partition 
number, M, and drawn in a log-log scale. In the top picture (A) 
the partition is performed in <f>— dimension only. The results of 
the partition in r\— dimension are given in picture B. The two- 
dimensional partition (n and <j>) is illustrated in the bottom part 
(C). We notice that the interacting system is obviously inter- 
mittent in one- as well as in two-dimensions. 



FM are calculated for the Pb+Pb events with the highest particle multiplic- 
ity (> 1200). The events with lower multiplicities are sampled, incompletely. 
The smallest one has a multiplicity of 410. As a result of this selection process, 
events at a rate ~ 0.002 per incoming primary beam are chosen for FM. Based 
on parameterization of charge-changing cross-section for heavy-ion interactions, 
we obtain apb+pb ~ 6.71 ± 0.21 barn p8[. According to the distribution shape 
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in the forward cone, the centrality of analyzed events can be estimated. The 
central ones contain two or fewer projectile fragmentations and have the high- 
est multiplicity. The other events are classified as peripheral or semicentral. 
Furthermore, the central events are characterized by total break-up into singly 
charged particles |^9| and supposed to be more favorable to produce QGP, as a 
result of the more nuclear matter they include. Besides the definition of interac- 
tion centrality, the multiplicity restriction can be employed to avoid any possible 
mechanical correlations Jl(| and to determine the kind of FM [jl2j . From these 
events, only the particles emitted within a predefined 77-interval, 2 < 77 < 6, are 
taken into account. This interval obviously contains the region of central ra- 
pidity. Due to the chiral symmetry breaking, the produced particles are mainly 
pions. Therefore, the restriction on rj is significant to study the deconfinement 
phase transition, which, in turn, is restricted to the produced particles. Here 
we employ an additional restriction on the considered phase-space. Its influence 
to get flat FM distribution has been discussed elsewhere jlO| . On the azimuthal 
space, there is no restriction, i.e. for each A77 group, <5</>-bins are allowed to take 
any real value within the available spectrum {0, 27r}. 

From the analyzed 380 events, 84 central events are selected according the 
criteria given above. Using the calculated cross section, this sample represents 
~ 22 ± 4% of all scanned events. The average multiplicity of their charged 
particles is ~ 1350 ± 50. From the multiplicity distribution per unit 77, we get 
average particle densityfJJ < dn/drj > ma x~ 540 ± 32. This results, according 
to the frequently used Bjorken formula p3| , energy density of 3.743 ± 0.324 
GeV/fm 3 . This value is evidently larger than the density required for the quark- 
gluon plasma formation |3C| . 

2.1.1 One-dimensional Analysis 

The pseudorapidity interval, 2 < 77 < 6, is successively divided into M equal 
bins. The multiplicity in each such bins is counted and then the corresponding 
g-order FM are calculated according to Eq. |l|. In Fig. ||A, FM are given in a 
log-log scale as functions of for the orders, q — {2, 3, 4, 5, 6}. The underscore 
in refers to the patitioned phase-space. The results of 77-partition are given 
in Fig. ||B. We notice that the dependence of FM on shows almost the same 
characteristics as in ^-dimension. We also notice that the experimental points 
can be fitted as straight lines. All lines have positive intermittency exponents, 
4>q, which evidently increase with increasing q. Generally, the slopes in 4>- are 
larger than the slopes in 77-dimension. 

2.1.2 Two-dimensional Analysis 

The multi-dimensional FM (</>, 77, • • •) |nj are suggested in order to study the 
interaction dynamics, to study the sources of multiplicity fluctuations, and con- 

1 The average multiplicity within m-th bin of size dr\ centered around the peak position, 
rjm, of dn/ffov-distribution. 
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sequently to clearly explain the power-law scaling. For the multi-dimensional 
partition, we should first utilize cither isotropical or anisotropical method. If 
A77, for example, is divided into the same number of bins as A</>, this method 
called self-affine partition (isotropical). Clearly, it leads to total partition num- 
ber, M 2 . In Fig. |^C, F q are drawn in dependence on M<j> ■ M v for the orders, 
q = {2,3,4,5,6}. In this picture, the 2D partition is performed isotropically. 
We notice that the experimental points can be fitted as straight lines. All lines 
have positive slopes, (j)^. These slopes increase with increasing q. Obviously, 
c/)^ are larger than the ID ones (<^, and </> v ). 

2.2 Intermittency Phenomenon 

To study the intermittent behavior observed in Fig. |^, we should first check 
whether this phenomenon can be understood by means of known physics, like 
conventional short-range correlations. Elsewhere we discussed other in- 
terpretations, like Bose-Einstein correlation (BEC), multi-particle cascade, ran- 
domization, hadronic Cerenkov-radiation, measurement bias, etc. It is known 
that the intermittent behavior is inconsistent with a-, Lund hadronization, ge- 
ometrical branching model, etc. 

2.2.1 BEC and Coulomb Reactions 

The contributions of BEC to the particle correlation have to be considered, espe- 
cially for identical particles. Furthermore, it is claimed that the intermittency 
is completely controlled by BEC and certain quantum statistical mechanism 
|1| |2j. Including BEC in the FRITIOF Monte Carlo code, is was possible 
to simulate the like-charged two-particle correlations . In spite of these re- 
sults, one should notice that the Dalitz decays and 7-conversations dominate the 
correlations, especially the lower-order ones of unlike-charged particles. Using 
the emulsion technique, one has almost no chance to directly determine neither 
the charges nor the momenta of produced particles. These measurements are 
essentially required to estimate the BEC. Nevertheless, I have invented three 
different methods to estimate BEC in the emulsion . It has been found that 
BEC nearly dominate the interparticle correlations especially within very small 
phase-space intervals^, i.e. very small relative momentum, q 2 . The intermit- 
tent behavior in very small phase-space intervals has been discussed elsewhere 
[ [To[ p"2[ , where not only a linear upwards trend was observed but a steep expo- 
nential one. This exponential increase can partially be understood according to 

2 Within such small intervals, the differences between the 77- values of the selected particles 
are correspondingly small. The relation between rj— differences and the relative momenta, 
q 2 = \qi — 52 1 2 1 is given as follows : 

q 2 = M 2 - {imf 

= [?it + 92t+ 2 ?lt92tcos(<?ii -<h)\ + 
[Ml t + M 2 t + 2M lt M 2t cosh(r?i - 772)] 

M is the invariant mass. 
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BEC. The 2D partition method represents another source for this exponential 
upwards increasing |34|, [55|. For relatively small M, as that in Fig. |2|, both of 
BEC and 2D partition method play a neglected role in FM p6[ . 

The Coulomb reactions represent an additional effect on FM. These final 
state interactions take place, when the produced particles are closely emitted. 
Depending whether the particles are bosons or fermions, these reactions increase 
oder decrease the relative momenta, Q. Coulomb final state interactions have 
therefore, non-neglectable effects on BEC and then on FM, espetially for small 
Q. The reactions between the identical poins and between the fragments rep- 
resent about 10% for Q < 1 MeV/c j57|. For increasing Q, the reaction rate 
decreases, exponentially p2[ . ft can be neglectable for few MeV/c. Therefore, 
for small M, as in Fig. |2|, the contributions of Coulomb final state interaction 
can be neglected. 



2.2.2 Physics of Intermittent behavior 

Here, we introduce another physics for the intermittent behavior. We will ex- 
amine whether it is able to describe the experimental data. In such as way, 
we simultaneously estimate the responsibility of the critical transition on the 
intermittent behavior, as introduced in Sect. [l]. 

1. Self-similar processes [j], |3S) (e.g. QCD parton cascading). At the 
critical point of phase transition, the correlation length will be on the 
increase and under the scale transformation the fluctuations are expected 
to display self-similarity. If the gaussian approximation can be utilized to 
describe the particle production and from Eq. we get 

%L = ^L.( q -l). (5) 

As a particular case of this general description, we refer to the mono- 
fractal behavior (see Eq. || below). The implementation of Levy index /x 
H,00§ reads 



(6) 



2" -2' 



H has a continuous spectrum within the region of stability {0,2}. The 



index, fi, allows an estimation of the cascading rate |3£|]. Obviously, Eq 
cannot be not applied in the tails of these distributions. Whereas, Eq 
is more effective. The two boundaries of Levy index, drawn in Fig. ^D, 
are corresponding to the degree of fluctuations in particle production as 
follows: 

3 In different experiments, such as that in |43[ p4| , it has been found that the index, fi, can 
be even outside this stability region. 
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(a) u — 2, minimum fluctuations from self-similar branching processes. 
This is a suitable condition to apply the scaling rule M , 

HO <7) 

Therefore, from Eq. [5] the anomalous fractal ratios, d q /d,2, are equal 
to q/2. For these conditions, lZ q are corresponding to the multi- 
fractal processes. Later, we will realize that the last scaling rule is 
not able to describe the experimental results given in Fig. |A,B,C. 

(b) n = 0, maximum fluctuations. Meanwhile beneath the critical point 
and according to QCD, the correlation length rapidly grows, it di- 
verges suddenly according to the statistical quantum physics. There- 
fore, near the critical point, i.e. neath the total randomization, it 
is expected that {d q \lZ q } — > const, and (f> q oc q — 1. At this point, 
the characteristics of interacting system can be compared with that 
of mono-fractal one EJ|, fig ]. This should not lead to the conclusion 
that the mono-fractal behavior alone is sufficient for the phase tran- 
sition. In Fig. ||, we will distinguish between these quantities and 
ascertain the mono-fractal behavior. 

Between these two limits, the approximation, Eq. ||, is no longer valid and 
should be replaced by Eq. ^. 

2. Critical phase transitions p, fiq (e.g. quark- hadron phase transition) . 
If QGP indeed is to be produced, the interacting system is supposed to 
suffer from thermal phase transition during its space-time evolution. At 
the critical point, <p q oc q — 1, d q = d2, and therefore, 1Z q /1Z = 1. In addi- 
tion to the model introduced in Sect. [|, there are auxiliary argumentations 
about utilizing the mono-fractality as signatures of the deconfinement. As 
the correlation length diverges at T c and if there are no long-range corre- 
lations, then d q are expected to be equal. Therefore, d q can be used to 
indicate the phase transition as follows: since the hadronization is sup- 
posed to occur, only if the QGP-matter suffers from a phase transition, 
then the final hadron-matter is expected to be intermittent with constant 
d q . In the other case, if the hadronization is ordinary a result of cascading 
processes, d q are linearly depending on q. Therefore, at the phase transi- 
tion, the ratio of q-th intermittency slope to the second one is depending 
on the orders, Q — l, (Eq. ||), 

^ = (9-1). (8) 

4>2 represents the dimension of the fractal sets in which the observed in- 
termittent behavior occurs @. In Sect. [2X2] , we will get <p2 — 1 for van- 



ishing fi. Obviously, this value is an evidence for mono-fractality. These 
patterns are also expected during the second-order phase transition from 
QGP to hadron gas El. 
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Then, we conclude that the dependence of 4> q /<j)2 on q — 1 reflects essential 
information about the reaction dynamics and about the physics of intcrmittency 
patterns. Fig. [| visualizes the differences between the possible sources of inter- 
mittent phenomenon discussed above. Meanwhile Fig. is used to encourage 
the conclusion, posted in this letter. 

In next sections, the g-dependency of the quantities, 4> q /4>2, d q /d2, lZ q /1Z, 
and C g will be examined and it will be shown how the intermittency and the 
fractal structure of multiplicity fluctuations are used as signatures for the de- 
confinement phase transition Q| . 

2.3 Dependency on the orders of FM 




B a 

*2L — «b— 




Figure 3: The relations between q-order FM and the second- 
order ones are drawn in log-log charts. The top picture, 
A, depicts the analysis in (f>— dimension. The relations in 
r\— dimension are given in picture B. The bottom picture, C, 
illustrates the analysis in two-dimensions. In all pictures, the 
experimental points are linearly fitted. We notice that the slope 
ratios, <f> q /<j>2, increase with increasing q and they are in <f>— space 
larger than in r]— space. Obviously, the 2D slope ratios are larger 
than the ID ones. 



In Fig. ^J, we noticed that the intermittency exponents, 4> q , are independent 
on M{^|0}. Therefore, the ratios, <f>q/<h, can directly be deduced from the 
relation F q (M^ v ^y) vs. F 2 {M {vW ), 
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\ogF q (M {lim ) = %L log F 2 (M {vm ) + c q . (9) 

This power-law is also valid for the multi-dimensional partition method p7{ . 
In a log-log scale, such relation in ^-dimension is depicted in Fig. ||A for the 
orders, q — {3, 4, 5, 6}. Fig. ||B illustrates the results in ry-dimension. We notice 
that the experimental points can be fitted as straight lines. Also, the slopes 
clearly increase with increasing q. The results from 2D analysis are given in 
Fig. |^C. For the comparison between 4> q /4>2 in ID and 2D, one should first 
use anisotropical methods, in order to overcome the additional effects of bin 



superposition especially in ^-dimension |1Q , |3J, 35 



2.3.1 Dependency of 4> q /4>2 on q 
Relation (^|) can be read as a power-law, 

F q ex if* (10) 

In Fig. [| it is clear to recognize that the powers, (3 q = (t) q /(f>2, increase with 
increasing q. Also the ID (3 q in the rapidity dimension are smaller than that in 
the azimuthal space. The 2D j3 q are apparently greater than the ID ones. 

Using a specific form of Ginzburg-Landau model to simulate the deconfine- 
ment phase transition in heavy-ion collisions led to the conclusion that j3 q are 
independent on the temperature T < T c . This guides the use of (3 q as a signature 
for the quark-hadron phase transition f45|| . 

& =(?-!)" (11) 

The powers, ^, describe the aftereffects of the phase transition from chaotic 
to coherent state. In our case, the aftereffects are the dimension ratios or ba- 
sically the fluctuations in the final state of particle production. Generally, v 
are independent on the details of interacting system. Therefore, v can be used 
to characterize the behavior of the measurable quantities at T beneath T c (e.g. 
(/> g /(/>2, d q /d,2 1 1Z q /lZ, £ g , etc.). In this regard, v are not a set of critical expo- 
nents in a conventional sense. Furthermore, they have a universal relevance p5| . 

In Fig. |]A,B the relation (3 q vs. q — 1 are illustrated in r\- and (^-dimension, 
respectively. As discussed above, (3 q can directly be determined from the slopes 
of the relation logF g vs. logi 7 ^ (Eq. ||). We notice that the experimental results 
cannot be predicted by any of the lines given by Eq. | (v = 1.0 @) or Eq. [l] 
(y = 1.304 ]45|). Comparing our results with previous experiments of KLM- and 
EMUOl-collaboration |^8|, ^9), we notice that even though all results can not be 
fitted by any of these ^-values, Pb+Pb results are relatively closer to the first 
line. Although, the other line is obtained from the Ginzburg-Landau description 
of a second-order phase transition [jiB] , it ist unable to fit any experimental data. 
As given in Sect. |2.2.2| , v = 1.0 characterizes the critical phase transition. This 
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Figure 4: The dependences of the ratios, f) q = 4> q /(j>2 on the 
orders, q — 1, for ID are drawn in A and B for 4>— and 
r\— dimension, respectively, where the experimental results (solid 
points) cannot be predicted neither by Eq. W nor by Eq. 0. In A, 
the triangles represent O-AgBr and p+AgBr at 200 AGeV 
and the squares are for S+Au at 200 AGeV ^j. Pb+Pb data 
are lower that the line with v = 1.0, meanwhile the previous re- 
sults are over the line with v — 1.304. The 2D analysis is given 
in C. d q /d2 are also studied in dependence on q — 1 and drawn 
in D (solid circles). The solid lines represent the two boundaries 
of Levy stable law. The dash-line is the implementation of Eq. ^ 
for fj, — +0.1. The open triangles represent the ID results. 



value is also obtained from 2D Ising model with second-order phase transition^ 
p2| . Ordinarily, v = 1.0 refers to the mono-fractality. Experimentally, different 
values are obtained for v But non of them can be related to the phase 

transition. 

We can conclude that the ID results in either ?/- or 0-dimension do not show 
clear signature for QGP. On the account of this restriction, one narrows the 
QGP signatures either in the rapidity or in the location. This restriction ist 

4 On the one hand, the appearance of intermittent pat terns during a phase transition has 
been proven in 2D Ising model long time ago [ pp| , [22] , pij . The intermittent fluctuations are 
characterized by a constant anomalous fractal dimension. On the other hand, we should utilize 
such statistical model to study the quark-hadron phase transition, since in lattice QCD at 
final temperature, which practically is close relevant to QGP, there is no direct calculation for 
the phase transition. Its characteristics can be obtained by numerical and analytical studies 
for statistical spin models. 
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that as it is the reason for this breakdown. For that reason, the search for QGP 
signatures should be performed through the posstible rapidity partitions and 
simulateneously over the possible locations. In such a way, we hope to disinte- 
grate the fine scale of QGP. 

The 2D dependence of (3 q on q — 1 is given in Fig. ^C. We notice that the 
relation (^), which has been suggested in B and successfully utilized in |5lj , is 
not able to fit the experimental data. This means that the self-similar branch- 
ing processes and consequently the minimal fluctuations are not responsible for 
the power-law scaling (Sect. |2.2.2| ). The argumentations given in Sect. E.2.l| 

are 

valid, if the intermittency is in fact the self-similar behavior. The dash-lines 
represent the power-fit according to Eq. |ll| with v = 1.304 JiH). The same pro- 
cess, but with v = 1.0, is represented by the dotted lines (Eq. |§|). We notice 
that the first power-fit is completely unable to describe the experimental data. 
Meanwhile, the second one gives a line very close to the experimental data (solid 
points). As given above, this result can be regarded as a signature of deconfine- 
ment phase transition. On the contrary to ID results, we get v = 1.0 for the 
2D dependency of (3 q on q. 

Furthermore, it is important to recognize that Eq. || and Eq. [ll] arc not 
belonging to Levy stable region. In the next section, we use Levy space to draw 
the q-dependency of d q jdi. In such a way, we try to reconfirm the previous 
results. 



2.3.2 Dependency of d q /d,2 on q 

From Eq. | and Eq. we get 

£ = ( g _l)"-i (12) 
Fig. ^D illustrates this relation. We notice that the ID resul ts (open tri- 



angles) are completely outside the Levy stable region (Sect. 2.2.2 ). In spite of 
this result, it is a significant finding to realize, that the values of dqjdi are ap- 
parently independent on q (d q /d2 f» 0.66 ± 0.04). Eq. || with /i = +0.1, is only 
able to describe the 2D results (solid circles). The values of 2D d q jdi are also 
constant (d q /d 2 « 1.0 ± 0.03). 

For small but distinctly non-zero the possibility of QGP formation mixed 
within cascading processes can not be withdrawn |3^| . If fi = 0, then 4> q /4>2 — 
q — 1 and correspondingly d q /d 2 = 1 (Eq. || and Eq. ^|). The predictions of 
Eq. U (dqjdi = 1), which, as given above, emphatically supports the conclusion 
of phase transition ||, are at the low boundary of Levy space. Once again, 
the obtained values of d q jdi are not consistent with being proportional to q as 
claimed in [pl| . On the contrary, they are equal. 
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Figure 5: In A fractal Renyi dimensions, !Z q /lZ, are de- 
picted versus q. For ID as well as for 2D data, we notice that 
IZq/lZ ~ 1. In B, the relations between log £ g and q in ID and 
2D are illustrated. In contrast to J5i| / we have here a clear ex- 
ponential dependence of ( q on q. In the bottom picture, C, £ 9 
are tested as function of the order, q, according to Eq. [zj. The 
linear fits are implied by the lines. 



2.3.3 Dependency of lZ q /1Z on q 

As discussed above, lZ q , one the one hand, measure the degree of randomization 
in the final state of particle production. On the other hand, they represent the 
fractal critical dimensions of the randomized intermittent clusters observed in 
the confined matter which is entirely bedded inside the deconfined one. As given 
in Sect. |l|, the thermal phase transition can be marked by constant 7Z q . The q- 
dependency of lZ q in ID and 2D is depicted in Fig. ||A. As expected, we notice 
that lZ q /lZ w 1. Near the critical point, d q /d,2 — > 1. Then the corresponding 
Renyi dimensions, 1Z Q — > 1Z, and 4> q /4>2 oc q — 1. 



As given in Sect. 2.3.4 below, near the critical point, the normalized expo- 
nents, £ g , show an abnormal behavior. They have an exponential g-dependency 
rather than a linear one. This will be utilized to distinguish between the scaling 



rules given in Sect. 2.2.2 
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2.3.4 Dependency of Q on q 

The scaling rule, Eq. |?|, is only valid, if the correlated g-tuplets are built up from 
single correlated particle-pairs and q — 2 non-correlated ones. Starting from the 
assumption given in |5l| that FM are defined as the integration of multi-particle 
correlation functions together with the lower-order ones which strongly effect 
the higher-order ones and defining the normalized exponents as 

c «=^/(0' (13) 

it is assumed that ( q is linearly depending on q (as given in Fig. 6 of Ref. ]51j| ). 
In Fig. ||b, we check this relation for ID as well as for 2D FM. It is clear to 
notice that Q q exponentially increase with increasing q. This result supports, 
on the one hand, the assumption that Eq. [7] is not able to describe the ratios 
(t> q /(j)2- On the other hand, it supports the conclusion that d q are constant for 
all orders of FM, as we noticed in Fig. ||. 

Using the scaling-rule, Eq. ||, to define the normalized exponents, we get 

C g = ^/(g-l), (14) 

Eq. |l4| has been tested in Fig. [|C. We find that £ g are linearly depending 
on q. This guides to the conclusion that meanwhile Eq. fj] failed to describe 
the Pb+Pb data, Eq. || is obviously able. Otherwise, we notice that C, q in ID 
decrease with increasing g, whereas in 2D there is a positive increasing. 



3 Summary and Conclusions 

In this latter, the non-statistical fluctuations in Pb+Pb collisions at 158 AGeV 
are investigated. First, I would like to sum up the results obtained so far. The 
analysis of FM in ID and 2D shows that this interacting system is obviously 
intermittent. The values of ID [3 q are considerably lower than the results ob- 
tained by KLM- and EMUOl-collaboration. Eq. (?] suggested in M and utilized 
in faild to simulate the dependency of ID and 2D (3 q on q — 1. Meanwhile, 
Eq. g results a line in ID more close to Pb+Pb than to the other data, the re- 
sults of 2D (3 q can effectively be fitted by this equation. The behavior described 
by Eq. || is predicted, on the one hand, if the interacting system suffers from 
thermal phase transition during its space-time evolution. On the other hand, 
the disability of this equation to perfectly describe the ID ratios leads to the 
conclusion that the Pb+Pb data are not uniquely mono-fractal. This results 
should not disturb the main conclusion posted in this letter, since on the basis 
of the fine scale of QGP, one can expect, that its signatures might be locally 
oriented. Therefore, the individual analysis in rapidity- or azimuthal-space may 
not be favorable to detect QGP. The 2D analysis is therefore, more effective 
than the ID. In addition to these results, we get d q /d,2 = 1 in 2D investigation 
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of the anomalous fractal dimensions. This value confessedly supports the hy- 
pothesis of the deconfinement phase transition. Levy index, fj, = 0.1, obtained 
for this 2D analysis, is an additional confirmation of this transition. We noticed 
that the fractal Renyi dimensions are constant (lZ q = 1Z). Also, in Eq. [l3| the 
relation between the normalized exponents, £ g , and the orders of FM shows that 
£ q exponentially increase. This dependence assists the conclusion that the rela- 
tion (^) is not able to describe the experimental data and therefore, it should 
be replaced by Eq. |L as done in Eq. 14. 



Finally, we come to the conclusion that the data sample used in this letter 
is intermittent in ID as well as in 2D. The intcrmittency ratios can be given by 
using simple relations, as that given in Eq. [§. Furthermore, according to certain 
statistical models, the g-dependency of anomalous fractal and Renyi dimensions 
and the index fj, evidently support the deconfinement phase transition. More 
data are however, needed to confirm the results and interpretations reported in 
this letter. 
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